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Abstract 

In a recent article (P. Wochner et al, PNAS (2009)) x-ray scattering intensity correlations around 
a ring, in the speckle diffraction pattern of a colloidal glass, were shown to display a remarkable 
~ cos(n<y9) dependence on the angular coordinate (p around the ring, with integer index n depending 
on the magnitude of the scattering wavevector. With an analytical derivation that preserves full 
generality in the Fraunhofer diffraction limit, we clarify the relationship of this result to previous 
x-ray studies of bond-orientation order, and provide a sound basis to the statement that the angular 
intensity correlations deliver information on local bond arrangements in a disordered (or partially 
ordered) system. We present a detailed analysis of the angular cross-correlation function and show 
its applicability for studies of wide range of structural properties of disordered systems, from local 
structure to spatial correlations between distant structural elements. 

PACS numbers: 61.05.cp, 61.43.-j, 61.43.Dq, 61.43.Fs, 82.70.Dd 
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I. INTRODUCTION 



{F{^))^ = ^l F{^)d^ (2) 







In a recent experimentji by Wochner et al, partially coherent x-rays with a wavelength of 
0.154 nm were used to investigate a colloidal glass composed of PMMA spheres of 117 nm 
radius, with particle suspensions of concentration exceeding the glass formation value. The 
scattering pattern on a 2D detector has the concentric ring structure characteristic of dis- 
ordered systems, and the speckle appearance resulting from the partial coherence of the 
undulator x-rays at the European Synchrotron Radiation Facility (ESRF). The authors of 
Ref.- have introduced the 4-point angular cross-correlation function (CCF) Cq{A) defined 
as 

where I{q, (p) is the scattered intensity, q is the magnitude of the scattering vector q, cp is 
an angular coordinate around the diffraction ring of radius q, and 

denotes an angular average around the ring.^^ 

The remarkable observation by Wochner et al} is that, at least for some q values, Cq{A) 
is well approximated by a cosine function of an integer multiple of the angle A, i.e., Cg(A) ~ 
cos(nA); different values of n are observed for different q. In particular, the experiment on 
PMMA spheres, for q ~ 0.04 nm~^, showed a very clean cosine behavior with n = 5. The 
authors recovered similar behavior from a numerical simulation, assuming that icosahedral 
clusters are preferentially formed locally, and computing the scattering intensity and its 
correlations for a cubic lattice of 8 000 such clusters, with random rotational orientation. 

In the following, we derive analytical expressions for the Fourier series expansion of the 
correlation function in the < A < 27r interval, from which the cosine-like behavior of the 
angular correlation function is related to the arrangement and orientation of bond angles 
and interatomic distances in the system in a completely general way. We limit ourself to 
the Fraunhofer diffraction case here and leave the discussion of the Fresnel diffraction to a 
forthcoming publication. One interesting aspect of this phenomenon is that it is essentially 
two-dimensional in character; in a disordered three-dimensional (3D) system, it appears that 
among randomly oriented local bond arrangements, the largest effects are expected from 
local bond arrangements where n-fold symmetry axes are (at least approximately) lined up 



with the direction of the incident x-rays. This leads us on the one hand to investigate the 
relationship to pure two-dimensional (2D) systems: and in fact extremely marked effects, 
unveiled by previous x-ray studies of bond-orientational order in liquid crystals^"- (especially 
hexatic ones), are strongly related to the recent results of Wochner et al. Our aim is 
to provide a sound basis to the statement that the angular intensity correlations deliver 
information on local bond arrangements in a disordered (or partially ordered) system. In the 
dilute limit (where local entities containing specific bond angles are separated by distances 
much larger than the bond lengths) the angular correlation function can be explicitly related 
to a bond-orientational order parameter, which generalizes the order parameter proposed 
for hexatic liquid crystals^ by Bruinsma and Nelson^. 

In this first paper of a series of two, we give a general theoretical treatment of the problem 
of the x-ray cross-correlation analysis (XCCA) in a partially disordered system. In the second 
paper, we will present the results of various simulations that demonstrate the general findings 
presented here. This paper is organized in the following way. In the next section a Fourier 
series analysis of the intensity angular correlations is presented. In the third section a general 
theoretical treatment of the CCF is given and the expressions for the Fourier coefficients of 
the CCF's are derived for the case of a kinematical x-ray scattering. In the fourth section 
the contribution of different terms to the CCF is analyzed. A special treatment is given to 
dilute and close-packed systems. In the fifth section we consider correlations in 3D systems, 
when the effect of the Ewald sphere curvature becomes important. Conditions, at which the 
angular CCF shows cos (nA) behavior with odd n-numbers will be analysed. The paper is 
completed by the conclusions and outlook section. 



II. FOURIER SERIES ANALYSIS OF THE INTENSITY ANGULAR CORRELA- 
TIONS 

We generalize the CCF defined in Eq. ([1]) by introducing the intensity correlations at two 
different values of the momentum transfer vectora^i qi and q2 

„ ... _ (/(gi, v)nq2, V + A))^ - (/(gi, v))^ ■ (/(g2, ^))^ 
"'"^ ^" (/(gi,v^)),-(/(g2,^)), ' 

where the averaging over the angle ip is defined in Eq. ([2]). In the next section we will 
show, that the magnitudes of the scattering vectors qi and q2 are, in fact, the values of the 
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perpendicular components of the 3D scattering vectors qi and q2. One can readily see that 
the CCF ([3]) can be rewritten in a slightly different form 

which shows that we are dealing with the angular correlation function of a normalized 
deviation of the intensity on the diffraction ring. Let us define this quantity for each value 
of the momentum transfer vector qj as 

■^^IJshM^SI}^, (5) 

{I[qj,^))^ 

and notice that this function has obviously vanishing angular average. The measured corre- 
lation function (|3]) can therefore be written as 

C,„,,(A) = {Di{qi,v)Di{q2,V + ^))^- (6) 

In order to understand what periodicity or symmetry this function of A may display, let us 
now proceed to its expansion into Fourier series in the (0, 2tt) interval 



n=— oo 

Q,.. = ^ C,„,,(A)e-'"^rfA. (7b) 

Here C'^^ q^ is the n-th coefficient in the Fourier series expansion of Cq^^q.^{/S.). Substituting 
now the expression ([6]) into Eq. f l7bp and following the usual arguments for the Fourier 
transforms of convolutions we gefr^ 



CI,,, = Driqi)D^{q2), (8) 

where D^{qj) are the Fourier coefficients of a normalized deviation of the intensity. One can 
see, that in order to calculate the Fourier coefficients of Cg-^^g^i^A), one may first calculate 
those of Di{qj,ip), i.e., D]{qj) and then take a product according to Eq. ([8]). Note, that 
the definition (|5]) of the normalized deviation Dj{qj,(p) implies that D^{qj) = 0. Since 
scattered intensities are always real quantities, it is also easy to show that Dj"'{qj) = D^*{qj) 
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and, therefore, Cq_^"'g^ = Cg*g^. According to these symmetry conditions Eq. (17a|) can be 
represented in the following form 

oo oo 

C,„,,(A) = 2 ^ Re K,,,e-^] = 2 J] | C^^^^^ \ ■ cos(nA + 7„), (9) 

n=l n=l 

7n = arg(C;^^^gJ, 

where the summation is performed over the positive integer numbers n. 
In the particular case, when qi = q2 = q, Eqs. ([8]) and ([9]) reduce to 

oo 

Cg{A) = 2J2 Cg cos(nA), (10a) 

n=l 

C^=\D^{q)\', C^>0. (10b) 

The general analysis of CCF's presented in this section and particularly Eqs. fllOa[ llObp 
explain a single cosine behavior of CCF calculated from experimental data in Ref.-. Clearly, 
a strong single cosine dependence of Cq{A) can be observed only for those values of q, 
at which one of the Fourier coefficients significantly dominates over all others. In the 
following sections we will show how such coefficients can be related to the structure and 
symmetry of the system. 

It is to be noted that Eqs. f llOa^ llObp also imply that the Fourier analysis of the CCF ([T]) 



investigated by Wochner et al} does not really contain additional information with respect 
to the Fourier analysis of the dependence of the intensity. Examples^^"- of analysis of the 
periodicity in the angular dependence of the intensity can be found in studies of hexatic 
liquid crystal phases^^, performed with incoherent sources. 

In the next section, we present detailed derivations of the CCF, based on the kinematical 
x-ray scattering theory. 



III. GENERAL THEORETICAL TREATMENT OF THE CROSS- 
CORRELATION FUNCTION 

We start our discussion with a simple scattering geometry depicted in Fig. [1] A coherent 
x-ray beam scatters on the disordered sample and creates a speckle pattern on the detector 
in the far-field regime. As a general model system we assume a 3D sample consisting of 
identical 3D local structures (LS) of arbitrary shape, random orientation and position in 3D 
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space (Fig. [2]). Such a model includes a variety of systems, i.e., clusters or molecules in the 
gas phase, LS's formed in colloidal systems (similar to Ref.^), protein molecules, viruses or 
complex biological systems in solution. 

The coherent x-ray scattering amplitude A{q) from such a sample can be described in 
the first Born approximation (or kinematical scattering) as 

^(q) = / p{r)e'^-^dr, (11) 



where p(r) is a total electron density of the system^^. For disordered systems under consid- 
eration this electron density can be written in the following form 

N 

p{r) = Y,Pk{r-Rk), (12) 

k=l 

where pfc(r) is an electron density of the k-th LS at the position (see Fig. [2]) and the 
summation is performed over all LS's. Substituting Eq. (fT2|) into Eq. (fTT]) we obtain for 
the total scattered amplitude 

TV 

A(q) = 5^e^'i-^^Afc(q), (13) 

k=l 

where Ak{q) is the amplitude scattered by one LS, 

^*:(q) = J Pk{v)e'^ ''dv, (14) 

and the integration is performed over the volume of each LS. Eqs. (|T3l) and ([T3j) express a 
simple fact that under conditions of coherent illumination the total scattering amplitude for 
each value of the wavevector q is a coherent sum of the individual amplitudes from each LS 
modulated with the corresponding phase term exp(iqR^), depending on the position of 
each LS. 

Using Eq. ( JT3l) . we can write the intensity scattered at certain momentum transfer value 
q as 

N 



J(q) = e^'^-(^'=^-^'=i)A*^(q)A,,(q) 

ki,k2=l 



ki,k2=l 
N 

E / j PlM)pk.{^2)e'''^-^--^dr,dr^. (15) 

ki,k2=l 




Here, the following notation for the radius vectors connecting two particles 1 and 2 in two 
different clusters ki and k2 was used 

R-fcUi = ^-^2,^1 + r2i, (16) 

where Rfc2,fci = Rfc2 ~ R-fci is the radius vector connecting different local structures, and 
r2i = r2 — ri is the radius vector connecting subunits inside LS's (see Fig. |2]). 

We decompose now the scattering vector q = (q"*",?^) into two components: q-*" that 
is perpendicular, and that is parallel to the direction of the incident beam (see Fig. |3]). 
We define the perpendicular component of the scattering vector q-*- in polar coordinates as 
q-*- = (g-*-, (p). We also define the perpendicular R^^^^ = R^ + r^;^, and the z-components 
^klki — ^k2M + -^21 of the radius vectors introduced in Eq. (fT6|) (see Figs. [T]and[2]). Using 
these notations for the vectors we can rewrite Eq. f[T^ as 

N 

/(q)= e-^"'-'-^'-^ / /p:^(r^g^)p,2(r2\g^)e^'^"-^4"irfr^rf4. (17) 

fcl,fc2 = l 

Here we introduced a modified complex valued electron density function, defined as 

PkA^iA') = I Pt.{ri.z)e-^^'^dz. (18) 

We want to note here that our treatment is quite general and is valid for both cases of 
wide and small angle scattering. In the first case, the effect of the Ewald sphere curvature 
[see Fig. [3](b)], that manifests itself by the presence of the exponential factors e~*''^'^''2.'=i and 
g-fj^'z in Eqs. (fT7|) and (fT8|) . may become important. This effect could break the scattering 
symmetry of a diffraction pattern, characteristic for the scattering on a positive valued 
electron density (Friedel's law) and may reveal additional symmetries that can be still hidden 
in the small angle scattering case. As it will be demonstrated in our model simulations this 
wide angle scattering geometry may become important for a scattering on atomic systems 
with local interatomic distances of the order of a few Angstroms. In the small angle scattering 
geometry with the scattering angles 2a « 1 we have for the values of the scattering vectors: 
q ~ 2A;a(l-aV6+---)> - '2,ka{l-2a'^ /?,+...)] ~ 2ka'^{l-a^ /?,+...). It is well seen from 
these expressions that the q^ component of the momentum transfer vector is proportional to 
the square of the small scattering angle a. It means that, in this situation, the z-components 
of the momentum transfer vectors are much smaller than their perpendicular components. 
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i.e., « q and can be neglected. In this limit we have a simplified expression for the 
intensity f lT7|) that does not depend on the z-component of the scattering vector q^ . For a 
real valued electron density PkX'^i) the modified electron density function ( ITS]) reduces to a 
real valued projected electron density of a LS 

PkX^i) = j PkM.z)dz. (19) 

This case of a small angle scattering is typical for scattering on colloidal samples with a 
typical distances between colloidal particles of few hundred nanometers as in Ref.-. 

According to Eq. ([8]), the Fourier coefficients of the CCF are determined by the Fourier 
coefficients of the normalized deviation D'j{qj). Direct calculations (see Appendix A for 
details) give for D^{qj) 

D7{qj) = riqf, q-)/l'{qf. q]), 0, (20) 

where the Fourier coefficients of the intensity I"'{qf, qf) are 

N 

r{qt,q]) = {^r E e-'^I-^'=-'=iL^^,,^(g/, g|), (21a) 

ki,k2=l 

LlM(^t,q]) = I ldrid4KM^^])P>^^i^2,q-)Jn{qt\^^^^^^ (21b) 

Here Jnip) is the Bessel function of the first kind of integer order n, and 0rx2i is the az- 
imuthal angle of the perpendicular component of the radius vector R^^ defined in Eq. flTB]) 
(see Fig.E]). 

From the derived expressions we can draw the following important conclusions. According 
to Eqs. fl7a] [7b| [8|, [201 121a[ I21bl) . the initial four-point correlation function Cg-^^g^{A) can be 
represented by its Fourier series expansion, where each Fourier coefficient is defined by a 
product of two 2-point correlation functions of the form f l21al I21bl) . corresponding to two 
different momentum transfer vectors qi and q2. The magnitude of n-th coefficient is defined 
by the Fourier coefficients /"(gj'-,g|) in Eq. fl21al) . which depend through L^j^ i^^{qj-,qj) on 
the internal symmetry of LS's as well as on the medium range order of these LS's in the 
disordered system. We will discuss the structure of these Fourier coefficients for certain 
scattering geometries in more detail in the next sections. 
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IV. CCF DECOMPOSITION: LOCAL STRUCTURE AND INTERPARTICLE 
SPATIAL CORRELATIONS 



In this section we consider more closely the contribution of different terms in the expansion 
fl21ap to the Fourier coefficients C^^ We consider here a particular case of a 2D system in 
a small angle scattering geometry (2a << f), when we can neglect the 2;-components of the 
scattering vectors and g| . In this case, the modified electron density— 'pkX'^i) is defined 
by Eq. f|T9|) . The sum in the expression (12 lap for the Fourier coefficients of intensity I"'{qj) 
can be split into two parts: 



N 



N N 
.k\=k2=k k\^k2 



(22) 



fci,fe2 = l 

where the first sum corresponds to the terms with ki = k2 = k, and the last one to the 
terms with ki ^ k2. 

a) Dilute systems 

It can be shown (see Appendix B), that for dilute systems, when the average distance 
D between the clusters is much bigger than the size c? of a single cluster, the contribution 
of the second sum in Eq. ( 122|) can become much smaller than that of the first one. In this 
situation the main contribution to the Fourier coefficients of CCF's will be determined by 
the first sum in Eq. (122|) that we will consider in detail below. As soon as for the first term 
in Eq. f l22l) ki = k2 = k and, therefore, Rfca.fei = and R^^ fci = ^21 we have an especially 
simple expression for the integral L^^^k^^k^qj) in Eq. (]21bp 




Ll{q,) = / / rfridr2p:(ri)pfc(r2)J„(g,|r2i|)e-'"<^-. (23) 



If all LS's have the same internal structure but are oriented and located in space randomly, 
the phase 0r2i in the exponent of Eq. ( |23l) can be defined as 

where (pk is the rotation angle of the fc-th LS with respect to the fixed angular orientation 
0^2^ of the LS in the origin of the coordinate system. In this case, for each LS the integral 
( l23l) can be expressed in the following form 

Ll{q,) = e-<^'=L"(g,). (25) 



Here contribution of each LS k is determined by its rotation angle (pk in the phase and the 
integral L'^iqj) is the same for all LS's 

L"(g,)= / /rfrirfr2p*(ri)p(r2)J„(g,|r2i|)e-^"<^°2i. (26) 



According to the structure of the integral L"{qj) its value strongly depends on the sym- 
metry of a LS and determines selection rules for the values n of non-zero Fourier coefficients 
^qi 92 ■ These selection rules can be used for identification of the symmetry of clusters in 
diluted systems. For demonstration, we calculate in Appendix C the integral L'"'{qj) for 2D 
clusters with the different rotational symmetry (see Fig. Hj). For example, for the cluster 
with 5-fold symmetry [Fig. Hl^d)] only n = lOi, {i = 1, 2...) will give non-zero contribution to 
the Fourier coefficients of CCF's. Note, that the Fourier coefficient with n = 5 is forbidden 
in this scattering geometry. 

In the limit of dilute systems, neglecting the second term in Eq. ( l22l) we have for the 
Fourier coefficients C'^^^g^ of the CCF 

N N N 

CI,, oc r*(gi)r(g2) « E E - = L''*{qi)L^^{q2) E ^^'^^'"'^ (27) 

k=l k'=l k,k'=l 

where (f)k',k = 4>k' — 4>k- One can rewrite the sum in Eq. ( |27|) in the following form 

TV 

E e*'"'^"'" = N^e'""^). (28) 

k,k'=l 

Here, the average over all local structure orientations is defined as 

(e^"-^) = j p{4>)e'''^dcp (29) 

and 

N 

p{<j)) = E '^('^ - ^k'^k) (30) 

k,k'=l 

is the probability distribution of angular orientations. The average {e^"'^) is, in fact, a 
generalization, for n 7^ 6, of the bond orientational order parameter, introduced for hexatic 
liquid crystals^ by Bruinsma and Nelsort^. 

Now, we will consider two different limits for possible orientations of LS's in 2D plane. If 
all LS's have the same angular orientation, i.e., all (f)k',k = 0, then the probability distribution 
function p(0) reduces to a delta function p(0) = 5(0). In this case of completely oriented 
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system (e*"'''^) = 1. It means that non-zero values of the Fourier coefficients C"^ will 
be determined only by the values of the Fourier coefficients L"'{qj) with the scaling factor 
proportional to A^^. 

In another limiting case, when all orientations are uniformly distributed in 2D plane, 
p{(f>) = l/(27r), and the angular average 

(e*""^) = ^ f e*"'^^ = <o. (31) 



2vr , 

has nonzero value only at n = 0. As the Fourier coefficient with n = is not contributing 
to the CCF [see Eqs. ([9], llOap ]. for a dilute 2D system with random orientations of LS's all 
Fourier coefficients C^^ of the angular CCF will be equal to zero. It also means, that in 
this situation it is not possible to determine the symmetry of LS's from the analysis of the 
angular CCF. This is similar to the situation in a small angle x-ray scattering (SAXS), when 
there is no preferential orientation in the disordered system. 

In the case of partial ordering, angular orientations of LS's can be described, for example, 
by the Gaussian distribution. Such situation may be realized when a disordered system is 
in an external field (magnetic, electric, etc), which drives it towards a more ordered state. 
In this case the probability distribution is given by 

p{<P) = l/(av^)exp[-0V(2fx2)], (32) 

where a is the standard deviation. For this partially ordered state the orientational order 
parameter {e'^"-'^) is equal to 

(e-^)=exp(^-lnV). (33) 

In this case the number of Fourier coefficients in the CCF is limited. The strongest contri- 
bution to the Fourier coefficients CJ^^ g^ is given by the lowest values of n and is stronger for 
more ordered systems (that correspond to lower values of a). 

b) Close-packed systems 

In the case of a dense system, when the average distance D between clusters is of the 
order of the size of a single cluster, the second sum in Eq. fl22|) can not be neglected. It 
can significantly affect the spectrum of the angular CCF. Taking both terms of Eq. (l22l) into 
account, the Fourier coefficients of the angular CCF can be written as the following sum of 
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four terms 



N N N N N N N N 

ki=k2 k:j=k4 fci=fc2 kz^k^ kij^k2 ^3=^4 k\^k2 k^^ki 

Term 1 Term 2 Term 3 Term 4 

A schematic illustration of correlations corresponding to these four terms is shown in Fig. |5l 
The first term is a product of two 2-point correlation functions, each of which correlates a 
separate LS with itself and, therefore, contains the information only on the internal structure 
of LS's. This term does not depend on the density of a disordered system and was considered 
in detail in the previous subsection. The second and the third terms are the products of two 
2-point correlation functions, one of which depends only on the internal structure of a LS, 
and another one correlates different LS's, separated in space. The fourth term is a product 
of two 2-point correlation functions each of which defines correlations of spatially separated 
LS's. Clearly, the terms 5*2 , 5*3 and 5*4 depend on the density of a disordered system. In 
the following, we consider in detail each term of Eq. ( IMll . 

Taking into consideration the results obtained in the previous subsection, the contribution 
of 5*2 and in Eq. (IMl) can be written as 

N N 

(e™^)L-(gi) Ll^,,{q2) + (e"^"^) Yl L^.^) L'' (^2) , (35) 



SJ} + S'' = N 



kij^k2 ks^ki 

where L^{qj) and ita^^i) defined in Eqs. IBip . and the angular average (e*""^) = 

N 

(l/N) ^ e™'^'= is defined similar to Eq. f l28|) . According to its structure the non-zero con- 

k=l 

tributions of these terms are defined by the same selection rules {L"'{qj) 7^ 0) as for the 
term S*". In addition, the value of non-zero contributions will be modulated by the orienta- 
tional order parameter (e™*^) and the spatial correlations between different LS's defined by 

Ll,k2(^,)- 

Finally, for the fourth term in Eq. f l34|) we have 

N N 

^4 = E E Lk:,k2i^i)Ll,,S^2), (36) 

k\^k2 k3y^k4 

where L^^ k2(lj) defined in Eq. (IBip . This term is determined only by the spatial corre- 
lations between different LS's, and can play a significant role in the close-packed systems. 
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We support our discussion by the calculations of the normalized contributions of all four 
terms in Eq. f lM|) to the Fourier coefficients (for the case Qi = q2 = g)- In our simulations 
we consider a 2D disordered system consisting of pentagonal clusters [see Appendix C and 
Fig. 131(d)]. A high density system with D/d = 1.5 and a low density system with D/d = 150 
were considered in these calculations, with the size of a pentagonal cluster d = 440 nm. Both 
systems contain 121 pentagonal clusters, and were characterized by the same set of in-plane 
angular cluster orientations {(pk}- These angles were defined by the Gaussian distribution 
(!32|) . with the standard deviation cr = 0.2 ■ 27r/5 (see Fig. [6]^a)). This distribution of angles 
covers all possible orientations for a 5-fold pentagonal cluster. 

The results of the calculations for q = 0.037 nm~^ are shown in Fig. MJo, c) (here, for 
clarity, only the first 25 Fourier coefficients are shown). In the case of a dilute system 
[Fig. Et^b)], the contribution of the Term 1 is strongly dominating over the contributions of 
all other terms. It reveals 5-fold symmetry by the presence of the Fourier coefficient with 
n = 10. Higher orders {n = 20,30,...) are not present due to the choice of the q value. 
For this system the contribution from spatial correlations between different structures is 
negligible. In the opposite case of the close-packed system [Fig. [6]^c)], the contribution from 
spatial correlations (Term 4) dominates over the contribution from the local symmetry of 
individual clusters (Term 1). The fourth term significantly modifies the frequency spectra, 
in particular by adding the coefficients which are not related to the internal structure of 
clusters. In this case, the Fourier coefficients with n = 10 and n = 12 are dominant in 
the shown range of the spectrum, but only one of them, with n = 10, is related to the 
internal structure of clusters. 

Our results show, that for dilute disordered systems the main contribution to the cross- 
correlation function Cq-^,q^{A) is defined by the local symmetry of clusters. For a partially 
ordered system, one can extract this information by analyzing Fourier coefficients of the 
CCF. For a dense system, the spatial correlations between clusters can become dominant, 
and their contribution to the CCF can not be easily separated from the contribution defined 
by the internal structure of clusters forming the system. 
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V. CORRELATIONS IN 3D SYSTEMS. EWALD SPHERE CURVATURE EF- 
FECTS 



In our previous discussion of scattering on 2D systems, we have seen that only even 
Fourier coefficients of the CCF have non-zero values. Here we will show, that non-zero odd 
Fourier coefficients can be present due to scattering to high angles on 3D systems due to 
Ewald sphere curvature effects. In this case full expressions [Eqs. fl21a[ I21b|) ] containing 
z-components of the scattering vector g| have to be analyzed. 

To simplify our discussion, we will consider here a 3D system consisting of identical 3D 
clusters composed of identical point scatterers. The modified electron density f lTS]) of a 
cluster can be defined in the following form 

Mr\ q]) = f{q,) ^ '^(^^ " (37) 

i=l 

where f{qj) is a form-factor of a scatterer, and A^^ is a number of scatterers in the cluster. The 
coordinates (r^, Zi) define the position of the i-th. scatterer inside the cluster k. Performing 
the integration in Eq. (]21bp gives 



Ns _ 

Ll,,M,q]) = \f{q,)f Yl ^"'^'"-'UqflKtl)^'^^ "^^^'^ (38) 

l,m=l 

where the summation over index I is performed over the positions of scatterers in the cluster 
fci, and the summation over index m is performed over the positions of scatterers in the 
cluster k2- Substituting this expression into Eq. fl21al) we obtain 

N Ns _ 

/"(?/,?!) = (^)"i/(g.)r E m 

k\ ,k2=l l,m=l 

We note here that for ri 7^ the terms with ki = k2 and / = m are equal to zero. Taking 
now into account that the terms with interchanged indices, i.e. ki, ^2 and k2,ki, as well as 
/, m and m, /, differ from each other by a change of the sign of -Z'^'^^ and by an additional 
factor (—1)", which arises due to the change of the phase + TT, we have for 

even values of n in Eq. (15^ 

nqt,q]) = 2{zr\f{q,)f J] J] COS (gJZ.t.J J„(g/|RSj)e "^-"^ (40) 

l<A:i<7V l<l<Ns 
ki<k2<N l<m<Ns 
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and for odd values of n: 



— 1Yh(D I / 

r(g/,g|) = 2(^rM/(?.)l' E E "'^^'^^ • (41) 

l<fci<Af l<l<Ns 
ki<k2<N l<m<Ns 

From the performed analysis we can see that, due to the curvature of the Ewald sphere 
(non-zero g| component), we obtain non-zero odd Fourier components of CCF in scattering 
on a 3D system. These components become negligibly small at experimental conditions cor- 
responding to the flat Ewald sphere, considered in the previous section. A detailed discussion 
of differences between correlation analysis of 2D and 3D systems, based on simulations, will 
be given in the forthcoming paper. 



VI. CONCLUSIONS AND OUTLOOK 



The basic results of this paper, Eq. (12 lap and fl21bp . are characterized by the following 
structure. 

1. They break up as a sum over LS pairs. Two points belonging to two LS's of a pair 
define a phase factor through the angle of the projection of their connecting vector on the 
(x, y) plane. 

2. Additional oscillating factors come from the the Bessel functions of integer order 
depending on the projections on the {x, y) plane of the connecting vector and of the scat- 
tering vector qj-; and also from the effective density p7c(r^,g|). Note that, in the far-field 
diffraction limit adopted here, for scattering at small angles (small g^), odd m values are 
strongly suppressed in comparison to even ones by the trigonometric pre-factors. This is in 
disagreement with the strong m = 5 components observed in the experiment^, and will need 
an additional analysis, for example, in the near-field scattering geometry. 

3. Classes of LS pairs for which the oscillating factors systematically have the same 
sign give the largest contribution to the sum for a given n. This is the case for the ki = 
^2 pairs, and the (purely two-dimensional) examples described in Appendix C show how 
the non-vanishing values of n are related to the rotational symmetry of the LS's around 
a common axis aligned with the direction of incidence. However, each LS contributes a 
value multiplied by a phase factor related to its orientation with respect to a reference 
direction in the plane; it is then easy to see that if the ensemble of illuminated LS's has 
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a completely random orientation around the n-fold axis, the sum vanishes: in this case, 
indeed, LS's rotated by Ti/n with respect to a given direction are as probable as those 
lined up in that direction, and their respective contributions cancel in the total result. 
This is in agreement with the concept of bond orientational order^. If a non-zero Fourier 
coefficient is observed, it implies either a preferential alignment along a given direction, 
either by a specific physical reason or, alternatively, because the ensemble of probed LS's 
is small enough to display pronounced fiuctuations from the average uniform distribution 
of orientations; another interesting possibility, in view of the imminent availability of free- 
electron laser sources, could occur if the acquisition time is short enough to provide an 
"instantaneous" view, without effectively performing a time-average that necessarily restores 
the equal probability of all orientations. This may indeed be already the case in experiments 
involving very slow dynamics, as may be the case in^. 

4. In a three-dimensional fluid, the order parameter defined above is contributed only 
by those molecules for which the n-fold axis is, at least to some degree of approximation, 
aligned to the direction of incidence. This probably explains why the observed Fourier com- 
ponents, especially in the intensity, but also in the CCF, are weak when compared to the ex- 
tremely marked ones observed in hexatic liquid crystals, which are stacks of two-dimensional 
manifolds^i"-. It is tempting to speculate that the subset of LS's with an approximate line-up 
of a symmetry axis, in a three-dimensional system, is always " dilute" , in the sense that it 
is constituted by a small fraction of the total number of molecules or clusters. This would 
allow the application of results obtained in this paper for the dilute limit also to systems 
which are, in the three-dimensional sense, close-packed. In the companion paper, simula- 
tions are performed also with the purpose of establishing the extent of the deviation from 
perfect alignment of the symmetry axis which is compatible with an observable contribution 
to the CCF signal. It is important to bridge the gap between a two-dimensional theoret- 
ical interpretation that seems to arise naturally from the experimental geometry and the 
three-dimensional isotropy of ordinary samples. 

There are various directions that future experiments may explore; it would certainly be 
very interesting to monitor the CCF signal in a system in which a controllable experimental 
parameter (e.g. temperature, an electric or magnetic field) may provide a way to vary the 
degree of alignment of a symmetry axis; or in which the bond orientational order is well 



16 



characterized. 
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It is to be noted that the trivial angular dependence due to the linear polarization of the 
incoming synchrotron beam, important at large scattering angles, must be removed from the 
intensity. 

Similar cross-correlation function, which includes an additional step of averaging over diffraction 
patterns, was introduced by Saldin et al.—. 

Here we used the fact that all functions in the integrand are periodic functions with the period 

271. 

In Eq. (jlip we have also tacitly assumed an infinite illumination region that also means an 
infinitely small speckle size. In practice, due to finite size effects (a finite size of the coherent 
beam, a finite size of the coherent area in the partially coherent beam, or a finite size of a 
sample) the size of a speckle is finite and is of the order Aq ~ 1/d, where d is a typical length 
at the sample position. However, these finite size effects will not influence our further general 
treatment of the CCF's. 

In this section, we omit the superscript _L, assuming that all vectors are deflned in 2D plane. 
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Appendix A 



Here we calculate the Fourier coefficients 



r{q^.qn = — /(q)e-^"^dv' (Al) 

/TT 







of the intensity scattered at certain momentum transfer vector q, defined in Eq. (fT7|) . The 
scalar product q^ ■ R^j^fc^ in the exponent of ( |T7I) can be written as 

q^ ■ ni^l^ = ■ iR^^ij ,os(^ - <p^^,,J, (A2) 

where g-*", ip are the polar coordinates of the perpendicular component of the vector q-*- [see 
Fig. El^a)] and [R^^^^J, 0r^2i are the polar coordinates of the perpendicular components 
of the vector R^ffci (see Fig. [21). Substituting this expression in Eq. ( ITTl) and using the 
Jacob i- Anger expansion^ of the exponential functions in series of Bessel functions J„(p) of 
the first kind of integer order m 

oo 



m=— oo 



we can write 

N 



.2 X 



nq^,ql = E e-''^'-'^^'^^ J J plM,q^)p,,{ri,q^)dridri 

ifY. (A3) 



m=— oo 



Integration over ip in (IA3p gives 
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{dip/2-K) exp[i(m - n)ip\ = 6m,n, 







where 6m,n is the Kroneker symbol. Substitution the result of this integration into flASP 
finally gives 

N 



ki,k2=l 



• q^ = (^)" E e-^''"^'=-'=i / pl^{ri,q^)p,,{r^,q^)dridv. 



~2 ^ 



^(.^|RS,l)e (A4) 
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It is clear from the definition (lAip that 

fci,fc2 = l 

(A5) 

These results imply that the Fourier coefficients D'i{q) of the normalized deviation, defined 
in dS]), can be written as follows 

I if n = 0. 

Appendix B 

Here, we will prove, that for dilute systems the main contribution to the Fourier coef- 
ficients of intensity I^{qj) is given by the first term in the expansion fl22|) . In the case of 
a dilute disordered system, when typical distances between LS's are much larger than the 
size of a LS's itself, i.e., |R^-'^^^^| ^ |rjj|, we can use the following approximation for the 
values of the Bessel functions in the integrals ^^.^(gj, g|) in Eq. (121bp : J„(gj|R^^ fci I) = 
J„(gj|Rfe2,A:i + r2i|) — Jn{.1j\^k2M\)- Furthermore, for large values of the argument of 
the Bessel function Jn{p), p ^ ('^^/S — 1/8), one can use the asymptotic expansion^ 
J„(p) ~ y^cos { P ~ ^ — j) ■ Taking all this into account we, finally, get for the inte- 
gral L^-^_^k2 ^'i- ^21bp in the case of a 2D system 




Ll^kM = I I cirirfr2p:^(ri)p.,(r2)J„(g,|R^:,,J)e "'^^.^i (Bl) 

-^)e^"'"^-^^n;i^.., (B2) 



cos 



where 



drpkir) 



In deriving Eq. flB2p we also used an approximation 0j^x2i ~ (/)r that is valid for dilute 
systems and means that this phase does not depend on individial orientations of LS's, but 
is determined by their positions in the system. 

We analyze now this asymptotic behavior of the second sum in Eq. fl22p . According to 
Eq. (]B2p . for the given qj value the function L'^^^kiili) decays as 1/ ^/ qj\^k2,ki\ with the 
increase of the distance RA:2,fci between the LS's. At the same time, it can be noted that the 
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N N 

sum ■■■ ill Eq. ( I22l) contains iV(A^ — 1) terms comparing to terms in the sum Yl 

ki^ki ki=k2=k 

However, the presence in Eq. flB2p of the exponential factor e*"'^^'=2.'=i with random phases 
0Rfc2fc^5 corresponding to a big number of LS's present in dilute system, will additionally 
reduce the contribution of the second sum in Eq. ( I22l) . Therefore, for dilute systems the 
dominant contribution to the Fourier coefficients C" „ will be defined by the first sum in 
Eq. ([22]) with h = k2. 



Appendix C 

We consider here simple 2D structures (clusters) with distinct rotational m-fold symme- 
tries shown in Fig. |H We define the electron density of a cluster as a real- valued quantity 
in the following form 

Ns Ns 

Pin) = MQMr -r^) = J^ Mqj)Siri - r,)5(</., - </.,), (CI) 
1=1 1=1 

where Ng is a number of scatterers in the cluster, fi{qj) is a scattering factor of the i-th 
scatterer in the cluster, r = (r, 0), r and (j) are the polar coordinates of a scatterer in the 
cluster. Using the definition (ICll) and performing the integration in Eq. ( I26l) we get 

s,t=l 

For an arbitrary cluster with m-fold rotational symmetry shown in Fig. H] the following 
assumptions are valid: = a, i.e., all scatterers in the cluster are located on equal distances 
from its center, (pi = 2tt /m ■ [i — 1) , i = l...Ns, where m = A^s is a highest order of rotational 
symmetry in the cluster, and we also assume fi{qj) = f{qj)- Using these assumptions in 
Eq. (lC2p . we derive the expressions of L'^iqj) for each of the clusters shown in Fig. |H 

1. m = 2, Fig.m^a): 

L"(g.) = / ^l/feOn-^nlO) + '/n(2ag,)] if n mod 2 = ^^^^ 
^ 1 if n mod 2^0 

2. m = 3, Fig.m^b): 

3|/(gj)|'[-fn(0) + 2J„(y3ag,)] if n mod 12 = 
^"fe) = \ -6|/(gj)|V„,(V3agj) if n mod 6 = 0, n mod 12 ^ (C4) 

other n 
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3. m = 4, Fig.ll^c): 

' 4|/(g,)P{J„(0) + J„(2ag,) + 2J„(y2ag,)} if nmod 8 = 
Mfi(lj)\Vn{2aqj) - 2J„(v^ag,)] if nmod 4 = 0, 

n mod 8 7^ 
other n 



(C5) 



4. m = 5, Fig.ll^d): 



5|/fe)|'{^n(0) + 2[Jn{Aiaqj) + Jn(A2ag,)]} if n mod 20 = 
-10|/(gj)nj„(Aiagj) + J„(A2agj)] if n mod 10 = 0, 

n mod 20 ^ 
other n, 



(C6) 



where = ^i(5 - ^5), = -^1(5 + ^5). 

5. m = 6, Fig. m^e): 

' 6|/(g,)P{J„(0) + 2J„(ag,) + J„(2ag,) + 2J„(y3ag,)} if n mod 12 = 
6|/(gj)p[2 J„(agj) + J„(2agj) - 2J„(y3agj)] if mod 6 = 0, 

n mod 12 7^ 

other n 

(C7) 



Equations (]C3IIC7I) define the selection rules which determine the contributions to the 
n-th coefficient C^^ related only to the internal structure of clusters. For instance, 
Eq. (lC4p means that the contribution from the internal structure of the cluster shown 
in Fig. m^b) to the Fourier coefficients with n = 6,18,30,... is defined by the func- 
tion = -6\f{qj)\'^Jn{V3aqj); for the coefficients with n = 12,24,48,..., L'^(gj) = 
6|/(gi)|V„(V3agj); for n = coefficient = 3\f{qj)\'^[l + 2JoiV^aqj)], while other 
Fourier coefficients do not contain any information on the internal structure of this partic- 
ular cluster. 
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FIG. 1: (Color online) Geometry of the diffraction experiment. A coherent x-ray beam illuminates 
a disordered sample and produces a speckle diffraction pattern on a detector. The speckle features 
are defined by the finite size of the beam, or the finite size of the sample and its microscopic 
configuration. The direction of the incident beam is defined along the z axis of the coordinate 
system. 
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FIG. 2: (Color online) Disordered sample consisting of clusters randomly oriented and distributed 
in 3D space. The radius vector Rfca.fci = Rfc2 ~ Rfci connects the centers Rfc^ and of two 
different clusters ki and ^2, the vectors ri and r2 define the positions of particles inside each 
cluster with the origin of a local coordinate system in each cluster ki positioned at its center. 
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(a) 



FIG. 3: (Color online) (a) In the general case, the cross-correlation function can be defined as 
an angular average over two intensity rings with different magnitudes of the scattering vectors 
7^ iQ-t- T^^^ perpendicular components of the scattering vectors q;;*- and are defined in the 
plane {kx, ky) in the polar coordinate system as qj*- = [q^, ip) and q^ = (g^, ip + A), (b) Scattering 
geometry in the reciprocal space. Here k™ is the wavevector of the incident beam directed along 
the z axis, k^*"* and k2"* are the wavevectors of two scattered waves with the scattering angles 
2ai and 2a2- The scattering vectors qi = {(lj^,qf) and q2 = (q2";^2) decomposed into two 
components: q^ that is perpendicular and q? that is parallel to the direction of the incident beam. 
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(a) (b) (c) 

f m ^ 

FIG. 4: (Color online) 2D particles with simple geometrical shapes owning different rotational 
symmetries (rotational axes are perpendicular to the plane of the figure): (a) 2-fold, (b) 3- fold, (c) 
4-fold, (d) 5-fold and (e) 6-fold. 
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Term 1 Term 2 Term 3 Term 4 



FIG. 5: (Color online) Schematic illustration of different types of correlations contributing to the 
Fourier components C^^ ^^ of the angular CCF corresponding to four different terms in Eg. 1341 (see 
text). 
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FIG. 6: (Color online) Normalized contributions of different terms to the Fourier coefficients 
at q = 0.037 nm^^. A 2D disordered system consisting of pentagonal clusters was considered, 
(a) Gaussian distribution of the in-plane angular orientations of the pentagonal clusters (with the 
standard deviation a = 0.2 * 27r/5). The blue lines bound a central angle (j) = 27r/5. (b) The case 
of a dilute system {D/d = 150). (c) The case of a close-packed system {D/d = 1.5). 
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